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ABSTRACT 

The formulation and solution of inverse problems for the estimation of 
parameters which describe damping and other dynamic properties in distributed 
models for the vibration of flexible structures is considered. Motivated by a 
slewing beam experiment* the identification of a nonlinear velocity dependent 
term which models air drag damping in the Euler-Bernoulli equation is investi- 
gated, Galerkin techniques are used to generate finite dimensional approxima- 
tions, Convergence estimates and numerical results are given. The modeling 
of, and related inverse problems for the dynamics of a high pressure hose line 
feeding a gas thruster actuator at the tip of a cantilevered beam are then 
considered. Approximation and convergence are discussed and numerical results 
involving experimental data are presented. 
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L. Introduction 


The purpose of this note is to illustrate and explain some of the ideas underlying the use of 
parameter estimation techinques in investigating damping and other dynamic phenomena in several 
classes of distributed models for flexible structures. We do this in the context of two specific 
examples drawn from experimental structures. In the first example we present techniques and results 
that can be used to study nonlinear aspects of viscous damping (nonlinear air drag) in slewing 
maneuvers with flexible beam like structures. We shall describe some fundamental questions, 
present a model for which an estimation problem is of importance, and then show how this inverse 
problem can be approximated for computational purposes. We do this in a weak or variational setting 
and give convergence arguments to provide a theoretical foundation for the numerical schemes we 
have used. These convergence results are then followed by presentation of a numerical test example. 
Our methods have proved useful in current studies with experimental data, the results for which will 
be presented in detail elsewhere. 

In a second example, we outline techniques that have been useful in developing accurate models 
for the dynamic effects of a flexible gas hose/tip mass/thruster apparatus when it is attached to a 
flexible beam to provide an active control system (the so-called "RPL experiment"). Since detailed 
mathematical arguments for this project are given elsewhere [BGRW], we shall only outline the 
model, the approximation ideas, and present a summary of numerical results obtained when using our 
general approach with experimental data. 

2u. Nonlinear Damping in Slewing Maneuvers 

In a series of papers [HJ], [JH], [JHR], Juang and his co-workers describe experiments carried 
out to demonstrate the feasibility of actively controlling (stabilizing) vibrations of a beam during 
slewing maneuvers. Experiments were conducted with a 1 meter steel beam and with a 3.9 meter 
aluminum honeycomb solar panel cantilevered in a vertical plane and rotated in the horizontal plane. 
Each was attached to a torquing motor at the hub of rotation or "root" of the beam. In typical 
experiments, the beams were slewed 30° to 45° in 1.5 to 4.5 seconds. Strain gauges were located at 
the root and at .22 and .5 of the length i of the beam. An angle potentiometer (at the root) measured 
angular displacements during the slew. These measurements (strain, y xx , and angular displacement, 
0) were used as feedback to the motor which was then used to suppress the vibrations via an LQR 
theoretical formulation for the feedback control laws. 

A series of slewing maneuver experiments were carried out in a laboratory at NASA Langley 
Research Center with possible effects due to air damping present. These were then followed by 
repetition of the experiments in a vacuum chamber. Experiment and theory were in good agreement 
in the vacuum chamber experiments, but there were significant discrepancies between the theoretical 
model based simulations and the experimental data obtained in the laboratory setting. It is important 
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to understand the variation in responses of controlled flexible structures in such differing 
environments since many future large spacecraft studies will, by necessity, involve model 
extrapolation and design based on laboratory performance only without the benefit of vacuum 
chamber comparisons. 

It has been suggested that the inaccuracies in the Juang, et. al. simulations most likely resulted 
from two types of model error: (1) the absense of a viscous damping term to represent the nonlinear 
air damping and (2) dissipation at the fixed end of the cantilevered beams was not included in the 
model even though there was some loose "play" in the clamped or "built-in" end of the beams. This 
latter mechanism for energy absorption presumably should be modeled by some type of nonlinear 
boundary condition in place of the usual no displacement, no slope boundary conditions: y = y x = 0. 



Here we focus on a nonlinear air damping component of the model as suggested by Juang, et. 
al., and show that one might effectively compute this using the experimental data in a least squares 
setting. We use a formulation proposed by Juang and his co-workers; this model can be derived by 
combining first principle energy considerations with laboratory findings on nonlinear drag forces. If 
0(t) and y(t,x) denote respectively angular displacement (from some reference angle) at time t and 
bending deflection along the beam at time t and position x as depicted in Figure 2.1 above, the model 
including actuator dynamics has the form (see [JHR] for a more detailed discussion) 


dPj 

3t 2 


d 2 0 


dt 


d0 dy d9 dy 

— — + -r— ) + C. {x— — + tr— 

3 dt at 4 dt at 


a 4 y 


(2.1) p— + px — + c,f(x — + ■£) + M* — + 3T> + El — = 0 


ax’ 


_ T d 0 r* dy , r* aw ay. , r* , ao ay, , 

(2.2) I R — - + px — — dx + c,f (x— - — + -=(-) xdx + c.{x— + xdx 

B dt 2 J ° at 2 J o 3 dt at Jo 4 dt at 


r* ... 3 y 


d0 dy 


( l „ r.. de . a y, 


= X(t) 


d0 d 2 0 

(2.3) T(t) = k 0 e a (t) - k — - k 2 — 

at dt 


(2.4) y( t ,0) = |^(t,0) = 0, 


(.,<) = 4 M) = o. 

dx dx 
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Here e a represents the applied armature voltage (this term involves strain feedback in the control 
problem), the terms involving c 3 and c 4 represent viscous damping and the coefficients p and El are 
the usual beam material parameters, linear mass density and flexural stiffness, and I B is the moment 
lo px 2 dx about the axis of rotation. The damping function f is assumed to have the form f(v) = vlvl 
for v in some bounded region. Of course, one can assume without loss of generality based upon 
physical considerations that f becomes bounded as v becomes large, 
p A number of state space formulations of this model can be investigated in the context of 

identification and control problems including the following: 

t • d0 

| (i) Both 0 and 0 = — are treated as states, 

(ii) The time history for 0 is assumed known, but 0 is treated as a state (if 0(t) is known 
from experimental data, it does follow that 0(t) can be effectively obtained by 

\ differentiation of the data), 

(iii) Both 0 and 0 are measured reliably so that each can be treated as a known quantity. 

I For our discussions here of estimation of the damping coefficients, we assume that (iii) holds so 

that one can combine (2.1) and (2.2) to also eliminate 0 as an unknown in the model. We do this 
i before formulating a least squares problem for estimation of the damping. We also include a 

Kelvin-Voigt material damping term in the beam dynamics equation. 

We consider then the coupled system (actuator dynamics plus transverse vibrations of the beam) 
for 0< x < t , t > 0 

(2.5) py tt + px0 + c 3 f(x0 + y t ) + c 4 (x9 + y t } + El D 4 y + c D I D 4 y t = 0, 

(2.6) I B 0 +J o pxy t( dx + c 3 f(x0 + y t )xdx + c 4 J {x0 + y t )xdx = k Q e a -k 1 0-k 2 0, 
with appropriate boundary conditions for t > 0 


(2.7) 

y(t,0) = Dy(t,0) = 0, 


(2.8) 

EIDVM) + c D I D 2 y t (t,-6) = 0, 

El D 3 y(t,£) +c D I D 3 y t (a) = 0, 


and appropriate initial conditions. Here and below we shall use the notation D = 3/3x. Assuming 
that 0(t) and 0(t) are known for 0 < t < T, where T is the duration of the experiment, we may solve 
for 0 in equation (2.6), subsitute this into equation (2.5), and obtain a single nonlinear partial 
differential equation for the bending deflection y. Upon doing this we obtain 

(2.9) py tt + Ay tt ) + c 3 (t,x,y t ) + c 4 J 2 (t,x,y t ) + g(t,x) 

+ c 3 f(x0 + y t ) +c 4 {x0 + y t } + El D 4 y + c D ID 4 y t = 0 

where 
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(2.10) 

4y u ) - - 

px r« 

I B +k 2 J ° 

ps y u ds, 

(2.11) 

^(t,x,y t ) i 

P x 

W 

f f(s0 + y t ) s ds, 

(2.12) 

^ 2 (t,x,y t ) 

P x 

*B +k 2 

r l • 

•I {s0 + y } s ds, 
J o 1 

and 




(2.13) 

g(t,x)= - 

PX fk 

° 

e a (t) * kj 0(t)}. 


The boundary conditions (2.7), (2.8) involving no displacement, no slope at x = 0 and no moment, 
no shear at x = l , can be equivalently written as 

(2.14) y(t,0) = 0, Dy(t,0) = 0, D 2 y(t,*) = 0, D 3 y(t,i) = 0. 

Using this system, we may formulate an estimation problem for the parameter q= (qi.q 25 q 3 .q 4 ) 
= (El, c D I, c 3 ,c 4 ), given observations e ; , i = 1, ..., M, for the root strain D 2 y(ti,0). The problem 
becomes one of choosing from some admissible set of parameters QcR 4 (we treat here only the 
constant parameter case - variable coefficients can be readily treated with appropriate modifications of 
the arguments - see [BC1], [BC2], [BCR], [BR1] for details) a parameter q* which minimizes over 
Q the least squares criterion 

M 

J(q)= £lDy(t i ,0;q)-e i r 

i=l 

where y(-, • ; q) is the solution of (2.9), (2.14) for a given set of initial conditions y(0,x) = <b(x), 
y t (0,x) = 'F(x). Least squares problems of this type possess a number of interesting aspects (infinite 
dimensional states, theoretical and computational ill-posedness, etc.) which have been discussed by 
us and others in a number of previous publications. Here we focus on one particular question: state 
approximation techniques and the convergence arguments for approximating parameter estimates 
(these convergence arguments are also an important part of theoretical results which guarantee a type 
of inverse method stability - i.e. continuous dependence of estimates on the observations - see [B]). 
Before turning to these convergence arguments, we rewrite the system in a variational form 
(conservative form) with states v = y t , w = D 2 y. 

We first rewrite (2.9) as 

(2. 15) pv t + £(v t ) + q 3 Jj(t,x,v) + J 2 (t,x,v) + g 

+ q 3 f(x0 + v) + q 4 {x0+v} + qj D 2 w + q 2 D 4 v = 0 

w t = D 2 v 

v(0) = 4 # , w(0) = D 2 <b, 
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(2.16) 

(2.17) 


where we seek solutions v e H 4 nH 2 L ,we H| with 

H[ = {<pe H 2 ((U) 1 9(0) = D(p(0) = 0} 

hi = {<p 6 H 2 (o,-e)icp(-e) = D<p(.e) = oj. 

Then an equivalent weak or variational form is: Find (v,w) e H^x H° satisfying the initial condtions 

(2.17) and 

(2.18) < pv t + £(v t ) + q 3 J 1 (t,-,v) + q^t.-.v), <p > 

+ < g(t, ) + q 3 f(x0 + v) + q^x© + v}, (p > + < q t w + q 2 D 2 v, D 2 <p > = 0 for (p e 

(2.19) < w t - D^, \jr > = 0 for \\r e H°. 

We next make the following observations: for each of the terms involving q ^ in equation (2.18), 
there is an analogue involving q 3 , the difference being that the q 3 terms are nonlinear, the q^^ terms are 
linear. Since our primary emphasis here is in presenting convergence arguments for estimation of the 
nonlinear damping coefficients and since these arguments readily extend to the case involving the 
linear q 4 terms, we, for the sake of exposition, shall assume throughout the remainder of our 
discussions that q^j = 0. Thus the system of interest to us is 

(2.20) < pv t + £(v t ) + q 3 7 r 1 (t,-,v) + q 3 f(x0(t) + v) + g(t,-), <p > 

+ < qjW + qjD^, D 2 <p > = 0 for cp e , 

(2.21) < Wj-D^, \|/> = 0 for\|/eH°, 

(2.22) v(0) = *F, w(0) = D 2 <D. 

We approximate the system via the usual Galerkin procedures by choosing finite dimensional 
subspaces c and H N c H° satisfying certain approximation properties (as N — » °°) to be 
specified below. The approximate system for v N 6 H£, w N € H N is given by 

(2.23) < pv”+ £(v** ) + q 3 7 1 (t,-,v N ) +q 3 f(x0 + v N ) + g(t,-), (p N > 

+ < q t w N +q 2 D 2 v N , D 2 cp N ) = 0 for <p N e H£ , 

(2.24) < W N- D 2 v n , \j/ n > = 0 for \|/ N e H N , 

(2.25) v N (0) = P£ v(0), w N (0) = P N w(0) 

where P£, P N are the orthogonal projections (in the H° norm) of H° onto H£, H N , respectively. We 
then may define a sequence of approximating parameter estimation problems consisting of minimizing 
over Q the criterion 

(2.26) J N (q) = Xlw N (t,0;q)- ei l 2 

i=l 

where w N is the solution of (2.23) - (2.25) corresponding to q. 
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The procedure for arguing convergence of the approximate parameters q N (a minimizer of J N in 
(2.26)) is now well-documented in a number of our previous papers (see [B] for a summary). A 
crucial step entails the following: If {q N } is an arbitrary sequence in Q with q N — > q in Q, one must 
establish the convergence z N -» 0 , u N -» 0 where z N (t) = v N - P£ v, u N (t) =w N -Pw with v N , w N the 
solutions of (2.23)-(2.25) corresponding to q N and v,w the solutions of ( 2 . 20 ) - ( 2 . 22 ) 
corresponding to q. We outline the arguments to establish this convergence. 

Using (2.20) - (2.22) and (2.23) - (2.25) with q = q N , we find for all cp e H£ 

<pzN,(p> = < V N-v t + (I-P£)v t ,<p> 

= < (I - PN ) v t + tiy t ) - iXy"), cp > 

+ < q 3 J 1 (t,-, v) - q 3 Ji(t, •, v N ) + q 3 f(x 0 + v) - q^f(x 0 + v N ), (p > 

+ < qjW - q^w N + q 2 D 2 v - q^DV'*, D 2 9 > 

= < (I - P£ )v t + L(-zf ) + £((I - P£ )v t ), 9 > + (q 3 - q 3 ) < Ji(t,-,v), 9 > 

+ q^ < 7 i(t,-,v) - ^(t,-,^), 9 > + (q 3 - q 3 ) < f(x0 + v), 9 > 

+ q^< f(x0 + v) - f(x0 + v N ), 9 > - < q^(w N - Pw ) + q 2 D 2 (v N - P£ v), D 2 9 > 
+ < qjW - q^P N w + qjD^ - q^D 2 P£ v, D 2 9 > . 


Thus we find (suppressing some of the obvious function arguments) 


<p Z N+£(z?),9> = <d-P{?)v t +4(I-P£)v t ),9> 

+ (q 3 - q 3 ) < Ji(v), 9 > + q 3 < ?i(v) - 7 X ( p£ v), 9 > 

(2.27) + q 3 < 7 X { Pl v) - ^i(v N ), 9 > + (q 3 - q 3 ) f(x0 + v), 9 > 

+ q^ < f(x0 + v) - f(x0 + P£ v), 9 > + q^ < f(x0 + P£ v) - f(x0 + v N ), 9 > 
- < qN u N + q£D 2 z N , D 2 9 > + < qjW - q^ 1 P N w + q^^ - q^D 2 P£ v, D 2 9 > 

for all 9 e HjJ. Also, for y e H N we find 

< u^, y > = < wN - P N w t , y > = < w^ - w t + (I - P N ) w t , y > 

(2.28) = < DV* - D 2 v, y > + <(I - P N )w t , y > 

= < D 2 z N , y > + < D 2 ( P£ - I)v, y > + < (I - P N )w t , y >. 


To obtain convergence estimates we make particular choices for 9 and y in (2.27) and (2.28). Let 
9 = z N e H£ in (2.27) and y = u N e H N in (2.28) and define 0? = D 2 ( P£ - I)v, 0^ = (I - P N )w t . 


1 2 2 

Then from (2.28) we obtain, using the inequality ab < — a + eb where e > 0 can be arbitrarily 

4e 


chosen. 
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(2.29) 


Ii|u N l 2 S elD 2 z N l + J- l u N l 2 + el<l 2 + ±lu N l 2 + eie?l 2 + -L|u N l i 
2 dt 4e 1 4e 2 4e 

= elD 2 z N l +— i u N l 2 + el0^ l 2 + eie^ I 2 

4e 1 2 


Before deriving an analogous estimate from (2.27), we require some hypotheses on the nonlinear 
damping function f which will in turn yield estimates on L and J v Recall that we can expect f(v) to 
behave like vlvl = (sgn v) v 2 for v in some neighborhood (not necessarily small) of the origin, while it 
becomes bounded for v large. Thus, for our theoretical considerations here we make the reasonable 
assumptions that f satisfies a Lipschitz condition 


(2.30) If®) - f®) I < K I ^ 1 for % 2 real 


as well as satisfying a boundedness condition 


(2.31) If® ISM. 


From (2.30) we readily obtain the estimates (we shall need these to use in estimates from (2.27) 
with (p = z N ) 

< f(x0 + v) - f(x0 + p"v), z N > < J K I (I - p£) v(t,x) 1 1 z N (t,x) I dx 

(2.32) 0 

* K I (I - p£ )v(t) 1 1 z N (t) I <|l(I- P*v(t) I 2 I z N (t) I 2 , 

where we have used I- 1 to denote either the absolute value or the H°(0,£) norm (the interpretation 
being clear from the usage). Similarly, we find 

(2.33) < f(x0 + PN v) - f(x0 + v N ), z N > S K I v N - PN v 1 1 z N I = K I z N (t) I 2 . 


Furthemore, using the definition (2.1 1) of we easily find 


< JjO.-.v) - ^(t.-.P^v), z N > < J* j ™* - J* K I (I - P*)v(t,s) I s d s I z N (t,x) 
< p[|(I-Pl)v I 2 + I z N (t) I 2 ], 


(2.34) 


where the constant |i depends in an obvious way only on £,K,p, I B , k 2 . With similar calculations 
we find 


(2.35) < ^(t,-, P£ v ) - Jj(t, • v N ), z N > < 2\i I z N (t) I 2 . 

Using the definition (2.10) of L we next find 
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< £(zf ), Z N > = J J + k J 0 psz t (t,s) ds Z N ( t ,x) dx 


(2.36) 


= P T 7" J xzN (t,x)dx 
I 0 + k 2 2 dt [ J 0 


-p Id N. . 2 

— - -<x,z(t)>. 


„N X v N 


< 4(I-P L > t ),z > < A- l(I-P-)v t 


N. N 

, ) V I I Z I 


B "-2 


(2.37) 


< £ [ I (I - P^) v ( I 2 + I z N (t) I 2 ] 


where the constant \x depends only on p, I B , k 2 , i. 

Finally, we can now obtain a desired estimate from (2.27) with cp = z N . Using (2.32) - (2.37) 
we find 

T 2 

Id , N .2 P N2 1 N2lN2 

2 dF plz TTF <X ’ Z > ' 


I B + k2 


+ p[ l ( I - P L)v t |2 + | z N | 2 ] + |q 3 -q^| M | Z N | 

(2.38) + I I p.[ I (I - Pg )v I 2 + I z N l 2 J + | q ^ \2\x I z N I 2 + lq 3 - q^ I mVT" 

+ lq*l| [|(I-Pl)vI 2 + Iz N | 2 ] + |qj| K lz N | 2 

+ lq^llu N IID 2 z N !-<q> 2 z N ,D 2 z N > + I 0^ I IdV*! 

where 6$ = {qjW - q^P N w + q 2 D 2 v - q£D 2 P£ v). 

Noting that < x, z N > < -J Z /3 I z N I, we find that 


2 3 


(2.39) 


h +k 2 


,z N > 2 > * 4- lz N | 2 = ± 


! B+ k 2 3 


i + 3yp-e 


[l 2 3 

since I B = J px dx = p^ /3 We assume that Q is a compact set and that the admissible parameter 


set Q entails the constraints q 2 > v for some positive constant v. We then have 
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v I D 2 z N I 2 < < D 2 z n , D 2 z n >. 


Using these estimates in (2.38) we obtain 

1(1 q |rr^ I 2 p ^ ^ 1 fJ 7 

id?[ p ' > Js( 1+tl )ia-p^)v,i 2 

+ i+P + i + 3k(l + I + kH |z N,2 

+ I[m 2 + M 2 <] I q 3 - q, I 2 + — I u 21 1 2 + e I D 2 z N I 2 
z 4e 

-vID 2 z N I 2 + — 10^ I 2 + e I D 2 z N I 2 
4e 

or, assuming that I q 3 - q^l -> 0, 

( 2 - 4 °) lit PlZ * ' VkT <X,zN>2 * »Cl z N I 2 + At) + (2e - v) I D 2 z N I 2 

. 2 

k N.2 
+ — I u I 

4e 

where (F (t) is bounded and — > 0 as N — > «> under the usual assumptions (see [B], [BCR],[BR1]) on 
the approximation properties of H[J and H N (i.e. that P£ -» I, P N -> 1 in the desired topologies). Foi 
example, approximations based upon cubic splines in and linear splines in H° would suffice. 
Finally combining (2.29) and (2.40) we obtain 

(2.41) lit P ,Z |2 " + <x > zN>2 + luN|2 £ (3e - v) I D 2 z N I 2 

+ I z N I 2 + y 2 I u N I 2 + G N (t) 

where G N is bounded with G N (t) — > 0. Integrating this inequality and using the fact that z N (0) = 0, 
u N (0) = 0, we find 


(2.42) 


1 p 1 A) i 2 - < x, 2 N «) > 2 + i u N (t) i 2 | + (v . 3e) ( , D y< 

2 L B + *2 J J° 

< [y, I As) I 2 + y 2 I U N (S) l 2 J ds + [lG N (t) I + 1 At) |] dt . 


(s) I ds 


Finally, using (2.39) and choosing e such that v - 3e = 8 > 0 in this last inequality, we obtain 
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(2.43) 


2 


1 + 3k 2 /p i 


I z N (t) I 2 + 


+ 5 f I D z (s) I ds < A(N) + 

J n 


lu N («)l 2 j 

J o [y 1 1 z N (s) I 2 + \ 1 u N (s) I 2 j ds 


where A(N) — > 0. Thus, by the usual Gronwall arguments we have I z N (t) I 2 — > 0, I u N (t) I 2 -» 0, 
and |q I D 2 z N (s) I 2 ds — > 0, where z N (t) = v N (t) - v(t) and u N (t) = w N (t) - P N w(t). Under 
appropriate convergence properties for P£ and P N , we may then use the triangle inequality to obtain 

I v N (t) - v(t) I -> 0, 1 w N (t) - w(t) I — » 0, f I D 2 v N (s) - D 2 v(s) I 2 ds — » 0. 

j 0 


We note that the above arguments can be used to establish convergence at each t of the strain in 
the H°(0,£) norm. If we use the root strain D 2 y(t,0) in the least squares criterion, a stronger strain 
convergence (pointwise in the spatial coordinate) is needed to complete the theoretical development. 
Arguments in the spirit of those given above can be made to give strain convergence in the H 1 norm 
which, of course, yields root strain convergence. Arguments of this nature have been given for 
similar problems elsewhere [BCK], [BR2]; they involve some technical detail and we shall not 
pursue the development here. Instead we turn to a brief discussion of some computational aspects of 
these schemes. 

The approximation and estimation schemes discussed above can be readily used to develop 
computational algorithms for estimation of damping coefficients (including those for the nonlinear 
viscous damping terms). We have developed and tested numerically some software packages based 
on these ideas. We are currently using the packages with experimental data provided by J. Juang; 
these results will be reported elsewhere. We close this section with a brief summary of findings for 
one of the numerical test examples we have investigated. 

A test example was considered for the system (2.9) with l = p = El = 1.0 assumed known and 
c 4 = c D = 0. We sought to estimate the damping coefficient c 3 from simulated data for the strain. 
That is, we chose a particular function y(t,x) = .5t 2 (x 4 - 4x 3 + 6x 2 ) as the true solution of the system 
equation (2.9) with an appropriate forcing function added to the equation. The true parameter value 
c; = 2.0 was used in (2.9) along with choices of 0(t) and e a (t) which were qualitatively similar to the 
corresponding experimental time histories available to us. For example, we used 


1 / 7 

0(t) = 4l--e sin(co n pt + <{>) with p = ^ 1 ~ C » <t> = Arctan p/£ 


and £, co n chosen so that 0 has a maximum amplitude of 4.5 at t = 1.75. The values k 0 = 52.136, k t 

to 


> = 1672.96 and I B + k 2 = 32.95 were used in (2.9) along with 0(t) obtained by differentiating 0. 

For the approximating system we used modified cubic splines for both and H N ; the basis for 
H£ was the usual cubic B-spline basis modified to satisfy the boundary conditions of while the 
usual B-splines modified to satisfy the boundary conditions of H| were used to generate H N . In each 
case the subspaces H£, H n had dimension N + 1. In the fit criterion we used the root strain y xx (t i; 0) 
at 29 equally spaced observations in the time interval 0 < t < 7. For the initial guess c^ = 2.5, a 
residual of 1789 is obtained. For N = 4, the scheme produced a converged estimate c 1.996 with 
a corresponding residual of .0088. A number of other test examples were studied with equally 
satisfactory findings. 

i 

Hose Effects on the Dynamics of the RPL Structure 

The RPL structure is an experimental apparatus which was designed and constructed at the 
Charles Stark Draper Laboratory in Cambridge, Massachusetts with funding supplied by the United 
1 States Air Force Rocket Propulsion Laboratory (RPL). Its primary function is to serve as a test bed 

w for the purpose of investigating control algorithms and instrumentation (sensors, actuators, 

* processors, etc.) for the large angle slewing of spacecraft with flexible appendages. It was designed 

to specifically incorporate those features which make control design for large flexible spacecraft an 
especially difficult and challenging problem. In particular, this includes light damping, high 
flexibility, a large number of, and closely spaced natural modes of vibration, difficult to model and 
coupled structural and actuator dynamics and dissipation mechanisms, etc. 


flexible hoses 



11 


The structure itself consists of four aluminum beams, each of length -6=4 feet, width b=6 
inches and thickness h = .125 inches, cantilevered in a symmetrical fashion to a central hub which is 
mounted on an air bearing table. The air bearing table allows for the near frictionless rotation of the 
entire structure about the vertical axis. Control actuation is achieved via nitrogen cold gas thrusters 
mounted at the tips of two opposing appendages. The other two beams are passive with masses at 
their tips serving only to maintain the over all symmetry of the structure. Nitrogen gas is supplied to 
the thrusters from storage tanks mounted to the central hub through stainless steel, wire mesh 
wrapped, flexible, high pressure hoses. Electro-mechanical valves control the expulsion of the gas 
from the thruster nozzles. Each appendage is instrumented with a linear accelerometer at the tip. Data 
from the sensors is recorded and control input signals are generated using a MINC 11/23 
microcomputer. 

Effective control design depends heavily upon the availability of a high fidelity model for the 
plant. In the case of the RPL structure, it is immediately clear that a model involving partial 
differential equations would be of some use. For the transverse vibration of the passive beams, a 
distributed parameter model based upon the Euler-Bemoulli equation together with appropriate 
boundary conditions describing the coupled motion of the tip mass and the rigid body rotation of the 
central hub would be adequate. For the active appendages (i.e., those with the tip thrusters) on the 
other hand, a more sophisticated model which also captures the coupled dynamic effects (i.e. 
additional mass, stiffness and dissipation, torsional motion, etc.) due to the motion of the flexible 
thruster hoses is needed. 

Since the transverse vibration of each of the individual appendages is decoupled, for our 
investigation here, we consider the problem of modeling the hose effects on the transverse vibration 
of a single cantilevered (i.e. clamped - free) beam (see Figure 3.2). 



Figure 3.2 

We describe a model which was suggested by S. Gates of the Control and Flight Dynamics division 
of the Draper Laboratory wherein the hose is treated as a damped linear harmonic oscillator that is 
rigidly attached to the thruster assembly at the free end of the beam. More precisely, the hose is 
modeled as a proof mass which reacts against the tip or thruster mass via an elastic spring and a 
linear, viscous damper (see Figure 3.3) 
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Figure 3.3 


Letting u(t,x) denote the vertical displacement of the beam at time t and position x, 0 < x < Z , 

3u 

and assuming only small deformations ( i.e. Iu(t,x) I < < i, I - 5 — (t,x) I « 1 ) 

C7X 

we use the Euler-Bemoulli equation together with Kelvin-Voigt viscoelastic material damping to 
describe its transverse vibration . That is 


-,2 

a u 


d 4 du 


-A 

a u 


(3.1) p — (t,x) + c D I — — (t,x) + El — (t,x) = 0 0 < x < £, t>0 


at 


0X 


a x " 


where p is the linear mass density, c D and E are respectively the material coefficient of viscosity and 
modulus of elasticity and I = bh 3 /12 is the second moment or moment of inertia of the uniform, 
rectangular cross section of the beam. 

At the free end of the beam, corresponding to x = i , the motion of the tip mass m T and the 
hose mass m H are given by 

■\3 


m 


(3.2) 

and 

(3.3) 


a 2 u 

T at 2 


a' 3u a 3 u n 

(U)-c D I — 3 - (t, t) - EI — (t,-C) 

D a x 3 at a x 3 


= c H (&t) - $.(t, l)) + k H (y(t) - u(t,f)) + f(t), 


t > 0 


m H — Y (t) + C H ( -^ (t) - If + k H (y(0 - ^)) = 0, 

dt 


t >0 


respectively where y(t) is the vertical displacement of the hose mass at time t measured from the 
equilibrium position, f(t) is the thruster force at time t and c H and k H are respectively the hose 
damping and stiffness coefficients. Assuming that the rotational inertia due to the hose - thruster 
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assembly is negligible, we also have the zero moment condition 


(3.4) c D I — 4(t A) + El i) = 0, t > 0 

dx dx 

at the free end. At the clamped end, we have the usual geometric boundary conditions of zero 
displacement, 

(3.5) u(t,0) =0, t > 0 

and zero slope, 

(3.6) (t,0) = 0, t > 0. 
dx 

Assuming that the system is initially at rest, we have the temporal boundary conditions or initial 
conditions given by 

(3.7) u(0,x) = 0, |^(0,x) = °> 0 < x < i 

and 

(3.8) y(0) = 0, §(0) = 0. 

Our primary concern here is the inverse problem which is naturally associated with the 
mathematical model given by equations (3.1) - (3.8) above. The physical dimensions and mass 
properties of the beam and thruster assembly and the elastic properties of the material from which the 
beam is made are known. Also, the thruster output can be experimentally calibrated. Consequently 
the parameters i , p, E, m T and b and h (and therefore I) are known. The input function f is given as 
well. However, the coefficient of viscosity c D and the hose parameters m H , c H and k H must be 
determined via an identification procedure. Recalling that the structure is instrumented with a linear 
accelerometer at the tip of the beam, we formulate the following inverse or parameter identification 
problem. 

Given a known input f(t) and corresponding measured output z(t) for t e [t 0 ,t 1 ], determine the 
parameters q = (c D , m H , c H , k H ) T in a closed and bounded subset Q of R 4 , which minimize the least 
squares performance index 
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d 2 

(3.9) J(q) * f I ~(t^;q)- z(t)| 2 dt 

J t. L 


at 


where u(-, • ; q) denotes the solution to the initial boundary value problem (3.1) - (3.8) 
corresponding to the choice of parameters q = (c D ,m H , c H , k H ) € Q. 

Implicit in the statement of the problem above is the well posedness of the initial boundary value 
problem, i.e. the existence, uniqueness and regularity of solutions to the system (3.1) - (3.8). The 
infinite dimensionality of the distributed state constraints necessitates the development and use of 
some form of finite dimensional approximation. Both of these issues are most efficiently addressed 
via an abstract, functional analytic formulation of the system (3.1) - (3.8). 

Define the Hilbert space H = R 2 x L 2 (0,i) endowed with the usual inner product and let V = 
f(C,ri,(P) e H: <p e H 2 (0,£), cp(0) = Dtp(O) = 0, T) =<p(£)}. Define the coercive bilinear forms c 
and k from V x V into R by c(cp,\j/) = c H (£ -(p(£))(X -\j/(£)) + c D I <D 2 tp, D 2 \j/> and k(<p,vp r ) = 
k H (C - <p(*))& - V(-t)) + El <D 2 tp, DV> for cp = (£, <p(i), (p) e V and <\f = (j X , V, 

the operator M e £(H,H) by M(£,ri,(p) = (m H £, m T ri, pep) and set F(t) = (0,f(t), 0 ) e H. Then the 
initial boundary value problem (3.1) - (3.8) can be rewritten in weak form as 

(3.10) < Mu u (t), cp > H + c( u t (t), $ ) + k(u(t), tp ) = < F(t), tp > H , t > 0, tp e V 

(3.11) u(0) = 0, u t (0) = 0, 

for u(t) = (y(t), u(t,£), u(t, -)) e V. 

Depending upon the degree of smoothness imposed upon the input f as a function of t (i.e. Lj, 
Holder continuity, H 1 , etc.), standard results from the theory of abstract parabolic equations (see [F], 
[S]) can be used to demonstrate the existence and uniqueness of solutions to (3.10) - (3.11) with 
varying degrees of regularity. 

We define finite dimensional approximations to the system (3.10), (3.1 1) using a cubic spline 

based Galerkin scheme. For each N = 1,2, ..., let {p^}^ 1 denote the usual cubic polynomial B- 

splines defined on the interval [0,-fc] with respect to the uniform mesh {0,-t/N, 2-E/N, ..., f j and 
which have been modified to satisfy Pj^(0) = DPj*(0) = 0, j=l,2,...,N+l. Set Pq = (1,0,0), Pj*= 

N N+l 

(0, Pj (£), p ) , j = 1,2,..., N + 1 and V = span {p. . The Galerkin equations corresponding to 

(3.10), (3.11) for u^t) e V N are given by 

(3.12) < Mu£ (t), Pj N > h + c (u^(t), pp + k(u N (t), pp 

= < F(t), > H , t > 0, j = 0,1,2,..., N+l 

(3.13) ^(0) = 0 upO) = 0. 
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For q = (c D , m H , c H , k H ) e Q, let u^t; q) = (y N (t;q),u N (t,£ ;q), u N (t,-; q)) denote the solution to the 
finite dimensional system (3.12), (3.13) corresponding to the choice of parameters qe Q and let 

0 2 u 

u (t;q) be defined analogously with respect to the system (3.10), (3.11). If (t,-6;q) in the least 

at 2 

a 2 u N 

squares performance index J given by (3.9) is replaced with (t, l ; q) and the resulting finite 

at 2 

dimensionally constrained estimation problem (henceforth referred to as the N" 1 approximating 
estimation problem) is solved, an approximation to the optimal parameters q is obtained. Indeed, it 
can be argued (using estimates in the spirit of those given in the previous section; see [BGRW] for 
the details) that if f is smooth and if {q N } N=1 is a sequence in Q with q N — » q° , 

a 2 N 2 

then * ; q N )-> u (• ; q°) in L (0,T; H) (or equivalently, that — (-,£;q N ) -> — ^ (•, i ; q°) 

at 3t 2 

in L 2 (0,T)) asN-)». This in turn can then be used to show that if q N is a solution (which can be 
shown to exist) to the N^ 1 approximating estimation problem, the sequence {q N }^ =1 c Q admits a 

_N: „ _Nj 

convergent subsequence {q } j=1 with q q as j -> and q is a solution to the original infinite 
dimensionally constrained identification problem. 



Figure 3.4 
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We tested our model and scheme using data produced by the following experiment carried out 
for us on the RPL structure by Dr. Michel A. Floyd, formerly of the Draper Laboratory. With the 
central hub held stationary and the structure initially at rest, a thruster line was pressurized to 300 psi 
and the corresponding jet fired for 50 milliseconds. Linear acceleration at the tip of the corresponding 
appendage was recorded over the period 0 to 5 seconds with a sampling interval of 5 milliseconds 
yielding 1000 measurements. With the accelerometer calibrated to 5 volts/g (g = 32 ft/sec 2 ), the 
observations produced are plotted in Figure 3.4 above. 

From an FFT of the data we found the first three frequencies to be approximately .75 Hz, 7.5Hz 
and MHz. A visual inspection of the data immediately reveals the modes at .75Hz and MHz. The 
mode at MHz is a torsional mode, or, more pecisely, a twisting of the beam about its longitudinal 
axis. The excitation of this mode results from a combination of factors including the presence of the 
hose and the fact that it is connected to the thruster assembly at the top rather than the center of the 
beam and the opening and closing of the thruster valves. Since the torsional vibration is at a much 
higher frequency than are either of the first two transverse vibrational modes, and since it is relatively 
rapidly damped, in our study here we simply treated it as noise and left it unmodeled. A more 
detailed discussion of the torsional effects and its coupling into the accelerometer measurements of the 
transverse motion of the beam can be found in [BGRW]. 

The beam is made of a grade of aluminum having linear mass density p = .027 slug/ft and 
modulus of elasticity E = 15.84 x 10 8 lb/(ft) 2 . We have l = 4ft and I = bh 3 /12 = 4.71 x 10 8 (ft) 4 . 
The mass of the thruster assembly m T was determined to be .149 slug. A hose pressure of 300 psi 
was determined to be equivalent to a force of .2971b. We have therefore 

r0. 2971b 0 < t < .05 

= *0.0 .05 < t < 5.0- 

The index of approximation N was taken to be 4 throughout. A detailed study of the convergence 
properties of our scheme using test examples and simulation data was carried out with the results 
having been reported in [BGRW]. 

We neglected the hose effects and internal damping (i.e., we took c D = m H = c H = k H = 0) and 
used the standard Euler-Bemoulli model for a cantilevered beam with tip mass to generate values for 
the linear acceleration at the tip. The resulting acceleration profile along with the experimentally 
observed data is plotted in Figure 3.5. If one is willing to accept the Euler-Bemoulli theory as a 
reasonable description for the transverse vibration of a long, slender, flexible beam, then it is clear 
from the figure that the hose effects are indeed significant and should be modeled. The residual was 
found to be 3.03. 
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Figure 3.5 

In applying our scheme we used the data over the interval 3.0 to 5.0 seconds (where the 
contribution from the torsional mode has been significantly damped) sampled at the rate of 10 
observations per second. By matching the first two modal frequencies of the data with the first two 
modal frequencies of the model we obtained crude estimates for m H and k H that could be used as 
start-up values for an iterative optimization procedure. Then, using our spline-based scheme along 
with an iterative Levenberg-Marquardt nonlinear least-squares routine to solve the approximating 
estimation problems to search over m H and k H , optimal values for the hose mass and hose stiffness 
were obtained. Taking these values along with c H = 0 as start-up values, we used our scheme again 
to obtain optimal values for m H , c H and k H . Repeating this general procedure, we obtained the 
optimal values for the parameters c D , m H , c H and k H given in Table 3.1. The corresponding tip 
acceleration profile is plotted in Figure 3.6. The residual was computed to be .70 - a clear 
improvement. 

c D (lb sec/ft 2 ) m H (slug) c H (lb. sec/ft) k H (lb/ft) 

127.40 .0801 .0078 .413 

Table 3.1 
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Figure 3.6 


In summary, our results demonstrate the feasibility and effectiveness of using distributed 
parameter systems to model structural dynamics. Our finite dimensional approximation methods 
provide a viable means for estimating unknown parameters appearing in the models which can not be 
explicitly measured experimentally. Finally, our general approach appears to compare favorably (see 
[BGRW]) with other identification procedures for the estimation of parameters in models for 
structural dynamics which are commonly used in engineering practice. 
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